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Abstract. We study THz-emission from a plasma driven by an incommensurate- 
frequency two-colour laser field. A semi-classical transient electron current model 
is derived from a fully quantum-mechanical description of the emission process 
in terms of sub-cycle field-ionization followed by continuum-continuum electron 
transitions. For the experiment, a GEP-locked laser and a near-degenerate op¬ 
tical parametric amplifier are used to produce two-colour pulses that consist of 
the fundamental and its near-half frequency. By choosing two incommensurate 
frequencies, the frequency of the CEP-stable THz-emission can be continuously 
tuned into the mid-IR range. This measured frequency dependence of the THz- 
emission is found to be consistent with the semi-classical transient electron current 
model, similar to the Brunei mechanism of harmonic generation. 
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1. Introduction 

Since the demonstration a decade earlier by Cook and Hochstrasser [1] of intense THz 
emission from air ionized with a two-colour (w -I- 2a;) laser field, this phenomenon 
remains in the focus of attention. This is motivated by its potential as a source of 
intense single-cycle THz pulses, enabled by the high generated THz-field strength and 
large bandwidth and facilitated by the absence of material damage. The THz-range 
covers a wealth of fundamental resonances (in molecules: vibrational and rotational 
resonances, in solids: phonon and plasmon resonances, impurity transitions), which 
opens a wide field of possibilities for fundamental material and device research [UI3] 
as well as possible sensor applications, including the identification of atmospheric 
pollutants and use in food quality-control, atmospheric and astrophysical remote 
sensing, and (medical) imaging with unique contrast mechanisms [4]. 

Another motivation is the the interesting underlying plasma dynamics leading 
to the THz emission. Both an a; — (2a; — ui) four-wave mixing mechanism based 
on stationary or nonstationary third-order susceptibilities [S], as well as a tunnel- 
ionization-induced micro-current mechanism [filH have been considered. According 
to the electron current model, the THz emission from plasma originates from sharp 
quasi-periodic changes of free electron density occurring at the peaks of optical field 
ionization in the tunneling regime. The induced electron current rapidly varies on the 
sub-cycle time scale and yields generation of many harmonics, also known as Brunei 
harmonics 011]. Similarly to HHG where the mixing of the laser driving field with 
its second harmonic allows breaking the symmetry between the two quantum paths 
and yields generation of even harmonics [lOj , the THz emission can be considered an 
even (0*^-order) harmonic. This particular mechanism of THz emission is thus closely 
related to two prominent phenomena. Above Threshold Ionization (ATI) and High- 
order Harmonic Generation (HHG), also based on quasi-periodic sub-cycle ionization 
followed by electron acceleration in the driving laser field. High-order harmonics are 
generated via transitions of the electron between a continuum state and the ground 
state, which can be efficiently controlled on the single-atom and sub-cycle level mm- 

Here we show that the low-frequency THz emission corresponds to transitions 
between continuum states and can also be controlled by shaping the laser field. From 
quantum-mechanical continuum-continuum transition matrix elements, we recover, 
using the Strong Field Approximation (SFA) with a stationary phase analysis, the 
semi-classical interpretation in terms of transient electron currents by Brunei 0 
and Kim et al. 0. Furthermore, we experimentally demonstrate control over the 
timing of sub-cycle optical field ionization relative to the sign and value of the vector 
potential of the pulse m- This is achieved by by fixing the carrier envelope phases 
(GEP) and detuning the frequencies of the two optical driving fields away from 
commensurability [14) . This results in a very wide continuous tunability of the THz 
frequency emitted by the laser-generated plasma. 

This frequency dependence of THz emission received little attention as mostly 
the commensurate a; -|- 2a; schemes have been used so far. Incommensurate frequency 
combinations have been realized by superimposing a broadband 20-fs pulse with its 
narrowband second harmonic; however the tunability range in this case is limited to 
the bandwidth of the fundamental |14j . Employing the relatively broad tuning range of 
the OPA, we obtained the initial evidence for the wide frequency-tunability of the THz- 
emission [ISlilS]. Also, a scheme similar to our experimental setup based on a near¬ 
degenerate optical parametric ampliher (OPA) superimposed with the fundamental 
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Figure 1. (a) Scheme of mixing two incommensurate frequencies for generating 
tunable THz emission, (b) Scheme of the atomic transitions in a strong laser field: 
transitions between a continuum and a ground state yield HHG, while transitions 
between continuum states yield low order harmonic generation (including THz). 
(c) Perturbative four-wave mixing scheme for THz generation. 


laser pulse has been implemented m- However, due to lack of CEP stabilisation of 
the driving laser, field-resolved detection of the generated THz pulse was impossible. 
Furthermore, in that work only the lowest THz-frequencies (< 3 THz) and only very 
small detunings from the exact degenerate case have been considered. 

Here, we demonstrate an incommensurate-two-colour scheme for THz generation, 
where the ratio a;i/u ;2 is not an integer, based on a nearly degenerate OPA pumped 
by a CEP-stable Ytterbium-based laser emitting at 1030-nm wavelength. Inspired 
by our initial experiments [mile], we present extensive experimental data, theory 
and calculations that agree with the measurement data. We analyse the spectral 
properties of the generated THz emission over a wide range of detuning frequencies 
from the exact commensurate case. With the CEP stabilised input pulses we are able 
to generated CEP-locked THz pulses and detect them in a field-resolved measurement, 
thus demonstrating their applicability for time-domain spectroscopy probing the 
amplitude and phase response of a material. Due to the long wavelengths of the 
constituent colours (1030 nm and 1800-2060 nm), our scheme is a promising route for 
a high-efhciency [IH], CEP-locked, broadband and widely tunable mid-IR/THz source. 

2. Quantum mechanical transient electron-current model for THz 
emission 

A proper quantum mechanical treatment of THz emission is desirable to support the 
commonly used semi-classical two-step model by Brunei [8] and Kim et al. [^. A 
very successful fully quantum mechanical model that describes emission of an atom 
exposed to a strong laser held is based on the SEA m, where the inhuence of the 
atomic potential on the electron in the continuum, as well as the influence of all bound 
states except the ground state are neglected. Recently full ab-initio quantum TDSE 
calculations were performed to calculate THz emission nnnu, but the drawback is 
that the physical nature of the underlying dynamics is difficult to interpret. Here 
we use an analytical SEA model to calculate the atomic emission and derive the 
expressions that are identical to those previously attained in a classical description 
by Brunei of recollision-free harmonic emission. 

An electron in a continuum interacting with the strong laser held can be described 
as a wavepacket, which propagates in phase space: 

^c(i) = -* [ H(i,<b)VL(ib)17o(ib,0)^oclt', 

Jo 


( 1 ) 
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where Uo{ti,tQ) and U{ti,to) are the field-free and full propagators between times 
to and ti and the operator VL{t) describes the interaction with the laser field. Here 
the velocity gauge is used. The interpretation of Eq. m, which has proven both 
very fruitful and accurate, is that tb are the instants of strong-held ionization and 
4’c{t) is a superposition of many wavepackets emitted into continuum at times tb- 
Calculations based on semi-classical wavepacket propagators or quantum-trajectories 
explicitly utilize this picture by representing quantum mechanical amplitudes as sums 
over contributions of wavepackets moving along classical trajectories [22]. Therefore, 
we express Eq. o as a sum over such wavepackets th), t), emitted by the parent 

neutral system at times th and propagated to t: 

^'c = • (2) 

tb 

Note that we formulate these equations for a linearly polarized laser held, so the 
electron dynamics are limited to a single dimension along that polarization direction. 

The emission of electromagnetic radiation is proportional to the dipole 
acceleration, i.e. the time-derivative of the dipole velocity dcc{t) = (Tc \p\ 4'c), where 
p = tv is the momentum operator (atomic units are used throughout the paper). With 
the dipole velocity of the whole wavepacket is given by the sum of dipole velocities 
for all birth times: dcc{t) = J2tbi tb 2 \p\ 4’tbi)- Expressing the wavefunction of the 
free electron born at instant % as a plane wave, iptb [’’(i; ^b), i\ = \/P {th,t,r)e'^^^*^’*’‘^\ 
where S = pr — E {t — t\f) is the action and p is the electron density, it can be shown 
(see Appendix), that the diagonal terms (fbi = tb 2 ) are the stationary points in phase 
space and thus dominate this sum. Wavepackets emitted at different times, fbi ^ tb 2 j 
are separated in space, i.e. have a reduced overlap and lead to rapid phase oscillations. 
We can thus simplify the dipole velocity to a single sum over birth instants: 

dcc(t) = ^{tftblpl'fptb), (3) 

th 

and express the continuum-continuum transition dipole velocity matrix elements as: 

(V'tb \P\ V'tb) = j drp (tb, t, r) V (r, t, tb) + 5d{t), (4) 


with the velocity v = dS/dr. Here, the first term describes the case when the 
wavepacket overlaps with itself and dd{t) provides corrections. These may arise if 
some wavepackets emitted at different times have met at some place in the phase 
space later. 

This is an important difference as compared to HHG, where the continuum 
and bound parts of the two wavepackets evolve differently and accumulate a phase 
difference that depends on frequency and ultimately leads to chirped emission of 
attosecond pulses. In C-C transitions, however, the requirement for the wavepackets 
to overlap yields the Brunei harmonic radiation which is unchirped. 

We now concentrate on the first term in Eq. (|4|) and rewrite: 



p{tb,t,r) dr 


Wit,tb) 


J p(th,t,r)v(r, tb,t) dr 
f p(tb,t,r) dr 

-V-' 

v(t,tb) 


'^W{t,t^)v{t,th) (5) 


as a sum over all ionization instants tb of the product of IE(t,tb), the norm of the 
wavepacket created at time tb, and r(t,tb), the average velocity of the wavepacket. 
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The wavepacket norm is time-independent and is given by the ionization 
probability at W{t,tyj) = r(tb)i5tb, where St^, is the saddle point region and T 
the instantaneous ionization rate. Thus we have: 


dccit) = '^T{ty,)S%v{t,t^) = 

th 


r(tb)u(t,tb)dtb- 


( 6 ) 


Here we see that dcc(t) is the driving term for the Brunei harmonics [8] and 
would correspond to the current density in the quasi-classical model of Kim et al. 
[5] if multiplied by the atomic density in the generating medium. The emitted 
electromagnetic field is proportional to the dipole acceleration 

dccit) =T{t)v{th = t;t) + J T ith)vith,t)dth. (7) 

The first term of this equation is proportional to the small initial velocity after 
tunneling: v (th = t; t) ~ 0, which is neglected in the usual three-step picture, and the 
second one is the force F [t, r(t, tb)] = v (th, t) acting at time t on the wavepacket at the 
position r(t,tb) (center of mass of the wavepacket). The force is given mainly by the 
laser field and the Coulomb potential of the ion F [t, r(t, th)] = Tlas (i) + ^Coui [?'(C ^b)]- 
In practice, the laser field usually dominates for electrons ionized in the strong-held 
regime, so the dipole acceleration is simply given by: 

dccit) = Fi^it) J r(tb)dtb- (8) 


The emission is simply proportional to the driving laser Held multiplied by the 
ionization steps f r(tb)dfb- The spectrum of the emission is obtained by Fourier 
transforming Eq. dH): 


dcci^) — Eiosi^j)E 


r ith) dtb 


(w), 


(9) 


which results in the convolution of the laser Held spectrum and the spectrum of the 
ionization steps. 

Now we apply the formalism derived above to the calculation of the emission 
spectrum in the case of a two-colour mixing scheme, which allows generation of even 
harmonics, including the THz sideband. Simulating our experimental conditions, the 
Held in the calculation is composed of a strong fundamental wave at the frequency 
Wp and a weaker signal wave at around half the frequency of the fundamental wave 
Eit) = -b where Epit) and E^it) are the envelopes of the pump 

and signal waves respectively. The direction in which the electron is emitted is given 
by the sign of the velocity Vdith) = —>■ oo,tb) at which the electron drifts after 

the laser pulse is over. Fig. [2] b) illustrates the asymmetric case of the two colour 
Held where all of the electrons ultimately fly in one direction. The build-up of this 
directional electron current is responsible for the emission of the low-frequency-THz 
waves (see Fig. [2]c,d). 

Tuning of the central frequency occurs when two incommensurate frequencies are 
mixed. This is illustrated in the right column of Fig. [2] for a 0.6w -I- w-combination, 
and compared to the commensurate case (O.Sw -I- uj) in the left column. While in the 
latter case, the growth of a directional net-current (Fig. [^Jd) leads to the emission of a 
low-THz-frequency Held (Fig. [^tjd), the modulation of the current direction resulting 
from the detuning of the signal-frequency by Aw leads to a 2Aw-shift of the center of 
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Figure 2. Simulation of the induced current, using Eqs. [6j{9] with 
r(£b) = exp[—2(2/p)^/^/(3|Fias(ib)|)] the classical expression = 

i for two-colour fields with two frequency ratios: the left and 

-b 

right-hand sides correspond to the commensurate and incommensurate frequency 
case, respectively. Shown are the considered driving electric fields (first row), the 
dipole velocity (Eq. 0, corresponding to the induced current (second row, blue 
line; the red line shows a low-pass filtered version), the low-pass filtered (retaining 
only the THz-sideband) dipole acceleration (Eq. 0, corresponding to the emitted 
THz-field (third row), and the spectral intensity of the emitted radiation (square 
of Eq. O (fourth row). The considered ionization potential is that of nitrogen 
molecules. Ip = 15.7 eV. The ojp-component of the driving field has a peak field 
strength of Fmax = 27 GV/m (corresponding to 1 x 10^^ W cm“^), and the 
cJs-component has Fmax = 13.5 GV/m and a relative phase of 7r/4. 


the emitted THz-spectrum (Fig. [2t,h). The shift of lowest order sideband frequency is 
then compared to the experimentally measured THz-emission frequency. Higher order 
sidebands are beyond the EO sampling detection range and were not measured in this 
experiment. 

This tunability is further analysed in Fig. Ob). The THz-emission appears as a 
zeroth order sideband alongside higher-order sidebands [5]. The THz sideband starts 
close to zero frequency at degeneracy (uig = O.hWp) and its frequency increases as the 
frequency of the signal wave is tuned. This effect can be interpreted as a temporal 
phase modulation, shifting the THz sideband towards higher frequencies as illustrated 
in the inset of the Fig. HI)- Although the described THz emission mechanism is 
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Emission frequency O/cCp 


Figure 3. Simulated power spectrum of the emission (square of Eq. |9]l as a 
function of the frequency cUg of the tunable signal wave, for the same conditions 
as in figure [2] 


a phenomenon of intrinsically non-perturbative nonlinear optics, the resultant THz 
frequency tunability follows the law = 2ws ~ (Fig- Eb), as is also the case in 
the perturbative four-wave mixing picture, illustrated in Fig [Tb) . We also predict 
the feasibility of generating weaker sidebands corresponding to 6, 8, etc. photon 
processes. This process clearly enables a tunability through the complete THz and 
mid-IR spectral regions |14) . 

3. Experimental implementation 

An Ytterbium-based laser amplifier and a near-degenerate OPA were used to generate 
two-colour pulses that were used for plasma excitation (Fig. |4|). The Yb:KGW 
regenerative amplifier is actively CEP stabilized [23] and is used to pump a two-stage 
near-degenerate OPA. First, part of the laser energy is split into three parts. One part 
is used to generate a broadband super-continuum by focusing it in a sapphire plate. 
The part of the so-generated spectrum around 680-720 nm is used as a seed for the 
first OPA stage. The second part of the beam is frequency-doubled and serves as a 
pump for the first OPA stage. The resulting difference-frequency (idler) wave from 
the first stage (1800 nm to 2100 nm) is then amplified as a signal wave in the second 
OPA stage, pumped by the third part of the fundamental laser beam. The resulting 
CEP of the OPA output is thus equal to the pump laser pulse CEP phase [23], up to 
a constant: ips = Tp + const.. 

At the output of the OPA, the pump wave at Ap = 1030 nm and the signal wave 
at As =1800-2100 nm are separated using dichroic mirrors and then re-combined after 
adjusting the divergence of the beams and the relative delay. The polarisations of each 
of the two colour-component pulses are linear and set parallel to each other. The phase 
delay between the two constituent colour waves is = :/5p(As — Ap)(27rc)“^ -I- const., 
and can thus be controlled via the fundamental laser CEP ipp. 
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Figure 4. Scheme of the experimental setup. THz emission was detected via 
electrooptic (EO) sampling setup based on GaSe crystal. The probe pulses 
were generated in an additional OPA and compressed down to 20 fs in a prism 
compressor. 


3.1. Observation of directional electron current 

We verified that the two-colour pulses used for THz generation induce a directional 
electron drift by measuring the photoelectron spectra produced via above threshold 
ionization (ATI) with a stereo-ATI spectrometer [25l [2^ . The electron spectrometer 
consists of a small xenon-filled gas cell (< 1 mbar) placed in a vacuum chamber with 
two identical time-of-flight (TOF) electron spectrometers oriented opposite each other 
in order to measure the photoelectron spectrum emitted in the opposite directions. 
The generated photo-electrons pass through small holes in the gass cell into the TOF 
arms. Each TOF arm consists of a field-free propagation tube (« 50 cm) and a 
micro-channel plate (MOP) electron detector. The linear polarisation of each of the 
two-colour driving wave is in the direction of the TOF detectors. The signals from 
the two TOF detectors are then acquired with multi-scaler card and processed using 
a computer. 

The ATI spectrum in one direction measured for different CEP phases of the 
laser pulse is shown in Fig. [5] a).This measurement is done with the OPA tuned at 
degeneracy (As = 2060 nm) so that the maximum directional drift is produced, as 
illustrated in the left column of figure [2] The asymmetry of the low-energy direct 
electrons is opposite to that of high-energy rescattered electrons (not considered in 
figure [2]) and partially reduces the asymmetry of the total electron yield. However, 
because the number of direct electrons is much higher than of the rescattered ones. 

















































CEP-stable Tunable THz-Emission from Controlled Plasma-Electron Bursts 


9 



Figure 5. Stereo-ATI electron spectrometry measurement results obtained with 
two commensurate-frequency pulses: Ap = 1.03/im and As = 2.06/i.m. (a) Electron 
energy spectrum (log scale) measured on one side of the stereo-ATI spectrometer 
as function of the fundamental laser CEP. (b) Total electron count in left and right 
directions of the stereo-ATI spectrometer obtained by integrating the electron 
spectrum of the left and right detectors. 


the total yield still exhibits a modulation depth higher than 15%, as shown Fig.[5]b). 
These findings prove that in our experimental conditions, the induced electron current 
is highly directional and we are working in the tunnelling regime, as assumed for the 
theoretical model. 

3.2. CEP-locked, tunable THz emission 

We now demonstrate the phase-locked temporal modulation of this induced electron 
current to achieve tunability of resulting CEP-locked THz emission. 

For THz-generation, the combined two-colour driving waveform is focused using 
an / = 7.5 cm parabolic mirror into ambient air. The field of the generated THz 
radiation is detected using electro-optic (EO) sampling. In general, the bandwidth 
/max in the EO sampling detection scheme depends on the probe pulse duration 
fmax ~ l/(2rprobe) and it is the main limitation in our setup. The probe pulse of 
Tprobe ~ 20 fs duration was generated by amplifying the spectral portion at around 900 
nm of the white-light super continuum in a separate OPA stage and then compressing 
the pulse in a prism compressor near transform limit (see figure |4]). The probe pulse 
was characterised using SHG Frequency Resolved Optical Gating (FROG) with the 
result shown in Fig. [51 

The phase of the THz wave according to the frequency dependence is v^thz = 
2(ps — ^p. Therefore, the degenerate scheme where fundamental and frequency-doubled 
pulses are mixed (with CEPs (p 2 ui = const.) provides a passive stabilization 

of the THz-pulse’s GEP, {(frHz = ‘2<Puj — ‘P 2 uj = const), similar to difference- 
frequency generation based on a second-order nonlinearity (y^) [21] • Since in our case, 
‘Ps = Pp + const., the CEP of the THz-pulse is the same as that of the fundamental 
laser pulse: pthz = :/5p-l-const. Because of the field-sensitive EO sampling of the THz 
emission, active stabilization of the laser CEP is thus mandatory. In other words, for 
non-CEP-stabilized driving pulses, the THz emission is still emitted, but their CEP 
is random from pulse to pulse. 

This is why, as shown in Fig. |7^), without active CEP stabilization of the laser 
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Figure 6. Temporal characterisation of EO sampling probe pulse using SHG 
FROG, a) measured FROG trace, b) reconstructed spectrum and group delay of 
the pulse c) measured temporal profile and calculated Fourier limit. 




Figure 7. Experimental THz-generation results with cjp +ajs two-colour driving 
fields {ljs = (cjp H- Au)/2), for various detunings created by combining 

the fundamental wave (Ap = 1.03 pm, 250 pJ) with the tunable signal wave 
(As = [2.06,1.96,1.84,1.80] pm, 20 pJ). (a) Electric field of the THz transient 
measured using EO sampling in case of locked CEP or free running laser, for 
zero detuning, Au = 0, i.e. exactly commensurate frequency ratio, (b) Temporal 
traces of the EO sampling signal for different detunings Au, (c) THz emission 
spectra of plasma calculated from the EO signal, represented on linear scale, (d) 
The same spectra as in (c), represented on logarithmic scale. Each spectrum and 
temporal trace was normalised to its maximum. The arrows in panel (d) show 
the theoretically predicted central frequency of the emission. 


pulse, the emitted THz waveforms for each laser pulse are averaged out in our EO 
sampling. By locking the laser CEP, we create reproducible THz fields that let us 
study the two-colour-driven plasma dynamics. 

The detected THz emission spectra as function of the signal wavelength As < 2Ap 
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are shown in Fig.[7]b). Despite the rather narrow bandwidth of our 180-fs pulses that 
were used to excite the plasma, inid-IR emission frequencies up to ~ 40 THz were 
measured as the result of detuning the OPA from degeneracy. The frequency of the 
THz-sideband shifts towards higher frequencies due to temporal current modulation 
when the OPA is detuned from degeneracy, as illustrated by the numerical simulations 
shown in figure [2] We can detect only the lowest THz-sideband in the spectrum due 
to the limited bandwidth of the THz detection apparatus. At higher frequencies, 
the probe pulse duration is significantly longer than the half-cycle of the field to be 
measured. The decreasing sensitivity of the THz detection setup at higher frequencies 
also contributes to the broadening of the spectrum for higher frequency detuning. 
The experimentally obtained frequency dependence agrees well with the semi-classical 
transient electron-current model simulations in figures [5] and [3J we observe a shift by 
Ri Ajz for a signal wavelenegth tuned to uig = (wp -I- Ai/)I2. 

4. Conclusion 

The tunability and CEP-stability of the THz emission from a laser-induced plasma- 
spark is demonstrated by mixing two incommensurate CEP-stable optical driving 
fields, forming a parametric waveform synthesizer producing shot-to-shot reproducible 
strong-held multi-colour waveforms. The measured frequency dependence of the THz 
emission is consistent with the transient electron current model we derived based on 
continuum-continuum transition matrix elements. 

Acknowledgements 

These studies were supported by the ERC project CyFi 280202 and Austrian Science 
Fund (FWF) via the SFB NextLite F4903-N23. S. Haessler acknowledges support by 
the EU-FP7-IEF MUSCULAR and the by the Austrian Science fund (FWF) through 
grant M1260-N16. 

Appendix: Stationary phase analysis 

An evolution of a continuum electron wavepacket can in general be described by 
the propagator CD. For convenience we will replace the integral by a sum of small 
wavepackets corresponding to different birth times t\^: 

ipti, (r, t, th) = -ill (t, th)V (tb)17o(ib, 0)4'oi5tb (10) 

where Sth is associated with the width of the saddle point region around th that 
satishes the stationary phase condition: 

r = ro-\- f v(t",t)dt". (11) 

•Itb 

Here tq is within the size of the initial state 4'o and is therefore negligible compared 
to the oscillation of the electron in the strong laser field described by the second 
term. This condition is valid beyond the SFA (55]. By coherently adding up the 
small wavepackets of different birth times we get the wavepacket of the electron 
in the continuum, Eq. ©■ The evolution of '0tb(^) is governed by both the 
ionic potential and the laser field. To calculate the emission, the dipole velocity 
dcc{t) = {'Sc \p\ 4'c) = Stbi.tb 2 \p\ V'tbi) has to be calculated. In case of a strong 
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laser field, the greatest contribution to the the momentum operator p is given by the 
derivative of the action S which is a rapidly oscillating function: 

nc 


Then the dipole velocity can be written as: 


v)= E 

thl ,th2 


^tbi 


dS (t,tb2) 


dr 


^tb2 


ibi ,ih2 


(13) 

(14) 


where v = dS/dr is the velocity of the electron at point r, pi = p{t,t\^i,r) 
and p 2 = p{t,th 2 ,'r) are the densities of the two wavepackets. The non-negligible 
contribution to the sum is given by the diagonal terms tbi = ib 2 that satisfy the 
stationary phase condition. By taking into account only these diagonal terms and 
calculating the acceleration of the dipole, we finally arrive at equation ([8]) which is the 
same as the one used in photocurrent calculations. 
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